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Abstract 

The quasiparticle concept is an important tool for the description of many-body systems. We 
study the quasiparticle properties for dilute Fermi systems with short-ranged, repulsive interactions 
using effective field theory. We calculate the proper self-energy contributions at order (fcp/A) 3 , 
where A is the short-distance scale that sets the size of the effective range parameters and kp the 
Fermi momentum. The quasiparticle energy, width, and effective mass to 0((/cf/A) 3 ) are derived 
from the calculated self-energy. 
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I. INTRODUCTION 



A relatively new approach to many-body phenomena is that of effective field theory 
(EFT). The EFT approach is designed to exploit the separation of scales in physical sys- 
tems |[], |3|. Only low-energy (or long-range) degrees of freedom are included explicitly, 
with the rest parametrized in terms of the most general local (contact) interactions. Using 
renormalization, the influence of high-energy states on low-energy observables is captured 
in a small number of low-energy constants. Thus, the EFT describes universal low-energy 
physics independent of detailed assumptions about the short-distance dynamics. Further- 
more, the application of EFT methods to many-body problems gives a consistent organiza- 
tion of many-body corrections, with reliable error estimates, and insight into the analytic 
structure of observables. EFT provides a model independent description of finite density 
observables in terms of parameters that can be fixed from scattering in the vacuum. In low- 
energy particle and nuclear physics, EFT methods have successfully been applied to mesonic 
systems ||], processes with one baryon ||, and more recently also to processes involving two 
and more baryons || [Tj. 

A description in terms of quasiparticles is central to many approaches to the many-body 
problem. It has proven useful in many physical systems from condensed matter physics to 
nuclear physics [f| 0. If the states of the noninteracting system are gradually transformed 
into the interacting states when the interactions are switched on adiabatically, the excita- 
tions of the interacting system near the Fermi surface (the quasiparticles) have the same 
character as the excitations of the noninteracting Fermi gas. This is the basic assumption 
underlying Landau's theory of Fermi liquids |10, llf . The assumption holds, e.g., for a dilute 



Fermi system with repulsive, short-ranged interactions. This system has been a benchmark 
problem for many-body techniques for a long time [12|, [13], |14|, I5| . In particular, the low 
density allows for a controlled expansion in kpa s where \zf is the Fermi momentum and a s is 
the S-wave scattering length of the fermions. In the framework of EFT, this corresponds to 
an expansion in kp/A where A is the breakdown scale of the EFT. For a natural system the 
scale A sets the scale of a s and all other effective range parameters. Initially, there were not 
many applications for this model system, but recently ultra-cold Fermi gases of atoms have 
been achieved |R| The dilute Fermi gas also serves as a simpler test case on the way 
to an EFT description of nuclear matter. The quasiparticle properties in the dilute Fermi 



gas to 0((kp I A) 2 ) were first calculated by Galitskii [18|. An incomplete numerical study 
of the quasiparticle properties to 0{ikpj A) 3 ) was carried out in Ref. [JTSJ . In Ref. [OT, the 
proper self-energy was calculated off-shell to 0((kp/A) 2 ) and the momentum dependence of 
the effective mass was studied in detail. 

In this paper, we perform the first complete calculation of the quasiparticle properties 
in the dilute Fermi gas to 0((kp/A) 3 ) in the low-density expansion. At this order both the 
S-wave effective range r s and the P-wave scattering length a p contribute. We calculate the 
quasiparticle energies e&, width j k , and effective mass M*. As a check, we also calculate 
the chemical potential fi and the energy per particle and compare with the known results 
T% O, O, EH)]. The outline of this paper is as follows. In Sec. II we repeat briefly the 



elements of an EFT approach to dilute Fermi systems. In Sec. [TTT| , we calculate the proper 
self-energy at order (kp/A) 3 . The quasiparticle properties to 0({k,F / A) 3 ) are derived from 
the self-energy in Sec. [IV[ In Sec. [V], we summarize and present our conclusions. 
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II. EFFECTIVE FIELD THEORY FOR DILUTE FERMI SYSTEMS 



In this section, we review briefly the EFT approach for dilute Fermi systems. We will 
restrict ourselves to repulsive interactions. An EFT for such a system is discussed in detail 
in Ref. [15|. In contrast to more traditional many-body approaches [12, 13, 14] the EFT for 



the dilute Fermi gas is fully perturbative. It is appropriate to consider a natural EFT for 
heavy, nonrelativistic fermions of mass M with spin-independent interactions whose strength 
is correctly determined by naive dimensional analysis. The scale of all effective range param- 
eters is set by 1/A, where A is the breakdown scale of the theory. For momenta C A, all 
interactions appear short-ranged and can be modeled by contact terms. Therefore, we con- 
sider a local Lagrangian for a nonrelativistic fermion field that is invariant under Galilean, 
parity, and time-reversal transformations: 



£ = ^ id t + —]^-^(i{j^) 2 + %;\(^(ipV 2 i{j)+ H.c. 

Z lb 



2M. 

+^vW-0W) + ... , (i) 

where V is the Galilean invariant derivative and H.c. denotes the Hermitian 

conjugate. Higher-order time derivatives are not included, since it is most convenient for 
our purposes to eliminate them in favor of spatial gradients. 

For convenience, we choose dimensional regularization with minimal subtraction as our 
renormalization scheme. Calculating the two-particle scattering amplitude in the vacuum to 
0((k/A) 3 ), the coefficients G%% can be determined by matching to the leading coefficients of 
the effective range expansion (see, e.g., Ref. ]n| for details). Thus, the Gn can be expressed 
completely in terms of the effective range parameters: 

= ~W' 2 = 0_ 2~' 2 = ~M ' ( ' 

where a s , a p and r s are the S- and P-wave scattering lengths and S-wave effective range, 
respectively. 

We now turn to the finite density system. First, consider a noninteracting Fermi gas. 
The single particle states of the system are characterized by a momentum k and a spin 
quantum number s. In particular, the ground-state of the system is the state in which all 
single-particle states with momentum less than the Fermi momentum kp are occupied and 
all other single-particle states are empty. Any excited state of the system can be created 
by removing particles with momentum less than kp and adding particles with momentum 
greater than fop. 

Now we turn on the interactions. We assume that there is a one-to-one correspondence 
between the states of the noninteracting Fermi gas and those of the interacting system if the 
interaction is switched on adiabatically. Note that this is not the case if the interaction is 
attractive. In this case, it is energetically favorable to form Cooper pairs and our calculation 
does not lead to the correct ground state. The elementary excitations of the interacting 
system correspond thus to the particle and hole excitations of the noninteracting Fermi 
gas, and are referred to as quasiparticles and quasiholes, respectively. Information about 
the properties of quasiparticles is contained in the proper self-energy £*(/co, k). In terms 
of diagrams, the proper self-energy represents one-particle irreducible and amputated^ self- 
energy insertions. The proper self-energy is related to the full Green's function G a p(k) by 
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Dyson's equation 0: 



G aP (k) = G%{k) + G aX Ck)^Ck)x,G^(~k) , (3) 

where k = (fc , k), a,/3, ... are spin indices, and the subscript denotes the noninteracting 
Green's function defined in Eq. (||) below. For spin-independent interactions this equation 
can be solved analytically to give: 

Gap(k) = ^ ^ . (4) 

PK fc -k 2 /(2M) -E*{k) 

Here, we only repeat the Feynman rules for the proper self-energy: The vertices that 
follow from the Lagrangian ([]]) are illustrated in Fig. [j]. For short-range interactions the 
direct and exchange contributions differ only by a spin degeneracy factor. As a consequence, 
it is convenient to use so-called "Hugenholtz" diagrams where the direct and exchange 
contributions are combined in local four-fermion vertices. The Feynman rules are 



1. Write down all Hugenholtz diagrams that scale with a given order in kp as determined 
by the power counting rules given below. 

2. Assign nonrelativistic four-momenta (frequency and three- momentum) to all lines and 
enforce four-momentum conservation at each vertex. 

3. For each vertex, include the corresponding expression from Fig. [I]. For spin- 
independent interactions, the two-body vertices have the structure {8\^5ps + 5\s5p~/), 
where A, (3 are the spin indices of the incoming lines and 7, S are the spin indices of 
the outgoing lines. 1 For each internal line include a factor iGo(k)\.y, where k = (ko, k) 
is the four-momentum assigned to the line, A and 7 are spin indices, and 

-r, ns ' r ( 0{k-k F ) 6{k F -k) \ 

= ^ U-kV(2M) + ^ + fc -kV(2M)-* e J ' (5) 

4. Perform the spin summations in the diagram. In every closed fermion loop, substitute 
a spin degeneracy factor — g for each 5\\. 






-iC 2 ^ -iC 2 k-k' 

FIG. 1: Feynman rules for the vertices arising from the Lagrangian (|l|). The relative momenta 
of the incoming and outgoing particles, are 2k and 2k', respectively. The spin indices have been 
suppressed. 



1 It is most convenient to take care of the antisymmetrization in the spin summations. This is discussed in 
rule 0. 
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5. Integrate over all independent momenta with a factor / <i 4 £;/(27r) 4 where d 4 k = dko d 3 k. 
If the spatial integrals are divergent, they are defined in D spatial dimensions and 
renormalized using minimal subtraction as discussed in Ref. [[OJ. For a tadpole line 
with four-momentum k, multiply by exp(ifco^) and take the limit rj — ► + after the 
contour integrals have been carried out. This procedure automatically takes into 
account that such lines must be hole lines. 

6. Multiply every Hugenholtz diagram by the appropriate symmetry factor. 2 

Furthermore, we repeat the power-counting rules that determine the power of kp with which 
a given diagram scales: 

• for every propagator a factor Mj k 2 F , 

• for every loop integration a factor of kp/M , 

• for every n-body vertex with 2i derivatives a factor fcp/(^ 2,+3n_5 ) , 

where A is the breakdown scale of the theory which is determined by the physics not included 
in the EFT (such as the mass of a heavy exchanged particle). Thus, kp has to be much 
smaller than A, which is just the case for low density systems with short-range interactions 
where p cx kp <C l/a 3 ~ A 3 . In the above power counting rules, the external momentum k in 
E* is also counted as order kp. This is appropriate, since the external momentum is converted 
to kp if observables like the energy density or quasiparticle properties are calculated. All 
proper self-energy diagrams have an overall factor of the Fermi energy kp / (2M) which we 
suppress in the labeling of different orders. The leading order contribution then starts at 
C(A;_f/A). As a consequence, the expansion of E* reads: 

E*(fc , k) = El + E*(k , k) + E* 3 (k , k) + . . . , (6) 

where the lower index indicates the corresponding power of kp/ K. In a homogeneous system 
E*(/c ,k) depends only on k = |k|. The expressions for the self-energy through 0((kp/A) 2 ) 
for ko = k 2 /(2M) (sometimes called "on shell") were first calculated by Galitskii ||18|| . The 



general expressions for arbitrary ko were obtained in Ref. 
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2 This symmetry factor is most easily obtained by resolving the Hugenholtz vertices in Fig. |l| into vertices 
with direct and exchange contributions. For a Hugenholtz diagram with n vertices this substitution 
generates 2 n new diagrams. Now identify the symmetry factor that corrects for the double counting of 
topologically equivalent diagrams among the 2™ diagrams. 
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(4) (5) (6) 



+ ..A ..o... ..o..... 

(7) (8) (9) 

FIG. 2: Hugenholtz diagrams contributing to S|(fcoj k). Solid lines indicate noninteracting fermion 
propagator @, while dotted lines indicate amputated fermion propagators. The expressions for 
the vertices are given in Fig. [l]. 

with v = k/kp. If the imaginary parts in Eqs. @ and (|10D are expanded around the Fermi 
surface at v = 1, the first nonvanishing contribution occurs at 0((1 — v) 2 ) as required by 
Luttinger's general theorem [^TJ. In the next section, we calculate the proper self-energy 
contribution £3 at 0((k F /A) 3 ) for k = k 2 /(2M). As we will show in Section fVJ, only the 
on-shell value of £3 is required for the quasiparticle properties. 



III. PROPER SELF-ENERGY AT 0{{k F /kf) 



In this section, we calculate the contribution to the proper self-energy at 0((kp/ A) 3 ): 
S3. Yt* 3 ^(ko, k) is given by the nine diagrams in shown Fig. 0. We will omit the calculational 
details in this section and show a sample calculation of diagram (7) in Appendix [A[ Diagrams 
(1) to (3) vanish when evaluated. Diagrams (4) and (5) can be related to the derivative of 
£2(^0, w ith respect to ko: 



ss (4) (fcb,fc) + ss (5) (fcd,fc) = -s; dj: f k 0,h) 



k °=lM 



(11) 



They do not contribute the quasiparticle properties at 0((kp/A) 3 ) as will be shown in the 
next section. Diagrams (8) and (9) can be evaluated analytically and give purely real results. 
For diagrams (6) and (7), we could not obtain analytical results; instead we have simplified 
the analytical expressions as far as possible and calculated them numerically using the Vegas 
Monte Carlo integration routine ||22|| . Below, we give parametrizations for the contributions 



of diagrams (6) and (7) to £3 from a fit to the numerical results. We have fitted the results 
for different external momenta k to polynomials in (kp — k) in an expansion around the Fermi 
surface. We have truncated the expansion at the order in (kp — k) at which adding another 
term does not improve the x 2 of the fit anymore. According to this criterion a polynomial of 
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seventh degree is sufficient for the real parts, while the imaginary parts require a polynomial 
of ninth degree in (/cp — k). The results for the contributions to the proper self-energy at 
0((k F /Af) are: 



ReE* (6) (|^) = 



I m ^3(6) k) 



I m ^3(8)(^V?' *0 
I m ^3(9)(^V?' *0 



- 1)(<? - 3)^(A; F a s ) 3 |o.5908 - 0.8829(1 - v) 
+0.2140(1 - vf + 3.7284(1 - vf - 2.7392(1 - vf 
-4.7509(1 - vf + 7.0977(1 - vf - 2.6709(1 - vf} , 

(g - l)^-(A; F a s ) 3 (o.7811 + 1.0589(1 - v) 
+0.4160(1 - vf + 1.9942(1 - vf - 3.1936(1 - vf 
-4.8352(1 - vf + 9.6964(1 - vf - 4.1846(1 - vf} , 



k\ {g + l) , u ,3/1 1 2 



(12) 



k"p 
2M 



k 2 F {9 - 1) ,3 2 fl 1 2 \ 



(13) 
(14) 
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2M 



7T Z 



ImS3 (7) (^g, fc) — 



+5.2913(1 - vf - 13.0528(1 - vf + 18.7234(1 - vf 
-15.4468(1 - v f + 6.7036(1 - v f - 1.1653(1 - v f^9(l - v) 

+ {-0.0700(1 - vf - 5.3866(1 - vf - 84.0177(1 - v f 
-721.4519(1 - v f - 3546.7045(1 - v f - 9902.6272(1 - v f 
-14617.8867(1 - vf - 8855.7811(1 - vf}o(v - 1) , 

S^^W [{-0.3835(1 - vf + 3.1281(1 - vf 
ZM 7r z LI 

-7.1868(1 - vf + 12.1392(1 - vf - 16.3114(1 - vf 
+14.4972(1 - vf - 7.2813(1 - vf + 1.5639(1 - vf}e(l - v) 

+ {o.3715(l - vf + 8.0530(1 - vf + 77.9691(1 - vf 
+532.4357(1 - vf + 2237.8794(1 - vf + 5462.8759(1 - vf 
+7100.2592(1 - vf + 3791.8271(1 - vf \o(v - 1)1 , 



0, 
0, 



(16) 



(17) 



(18) 
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where v = k/kp. In Eqs. (0), (|13|)> (0)) an d Q17D , we give the numerical constants to four 
decimal places to ensure that the parametrizations differ from the exact numerical values by 
less than 1%. Also, note that the parametrizations for the imaginary parts start with the 
(1 — v) 2 term and satisfy Luttinger's theorem plfl by construction. 



IV. QUASIPARTICLE PROPERTIES 

The elementary excitations of the interacting system correspond to the particle and hole 
excitations of the noninteracting Fermi gas, and are referred to as quasiparticles and quasi- 
holes, respectively. Information about the properties of quasiparticles is contained in the 
proper self-energy S*(/c , k). The singularities of the full Green's function determine the 
excitation energies e k of the system and their corresponding widths ^ k . In the quasiparticle 
approximation, the full Green's function ([|) can be written as 

17 

G(k ,k) = + regular terms , (19) 

k - e k - 27 fc 

where Z is the wave function renormalization. To determine e k and j k , we have to solve the 
equation: 

k 2 

tk + ilk- 2^ - S*(e fc + i7/ c ,/c) = 0, (20) 
while the wave function renormalization Z is given by 



1 <9E*(fc ,fc) 
dk n 



(21) 

ko=tk+ilk, k=k F 



The effective mass M* of the quasiparticle is defined using the group velocity of a quasi- 
particle state at the Fermi surface; that is the slope of the excitation energy at p = kp 



M* = k F 



f de k 
V dk 



(22) 



k=kp/ 

The effective mass can be extracted from the heat capacity of the system in the zero tem- 
perature limit 0: 

C v k 2 B TM*kp 

— = as T -> . (23) 

Various other observables can be computed from these results as well. The chemical po- 
tential fj, describes the minimal energy required to add an additional particle to a system. 
Consequently, it is given by the quasiparticle excitation energy at the Fermi surface: 

V=tk\ k=kF - (24) 

Another observable, which allows a simple check of the calculated excitation energies, is 
the energy per particle. One can compute the energy per particle directly from the proper 
self-energy £*(po,p)- It is easier, however, to use thermodynamical identities. The chemical 
potential at zero entropy S is related to the exact ground-state energy E by the thermody- 
namic identity: 

A*=fHl. at 5 = 0. (25) 



KdNJ 



v 



S 



Integrating this equation at a constant volume V and zero entropy gives: 



E = / dN'fji(3 = 0, V, N') . (26) 
Jo 

Since the N dependence of /i is hidden in the Fermi momentum k F through the relation 
k F = (<6-K 2 N/gV) 1 ^, the energy per particle is easily evaluated by multiplying each term in 
the chemical potential of order (kp) x by a factor of 3/(3 + A). 

To obtain all these quantities, we solve Eq. ( p0|) perturbatively in hp/ A. We therefore 
expand e k and 7^ in powers of kp/ A as well: 

Cfc = Cfc.O + e k,l + Cfc,2 + ■ ■ ■ , 

Ik = 1k,0 + 7fc,l + 1k,2 + ■■■ ■ (27) 

Inserting Eqs. (§) and (p7|) into Eq. (|2~0"D and requiring Eq. (|20|) to hold for every order in 
hp/ A, we obtain the matching equations. The solutions through order (kp/A) 2 were first 
obtained by Galitskii |18| (see also Ref. |J). For completeness, we quote their results: 

k 2 

e k,0 - J 7fc,0 - , 



Efe,i = S* , 7 fcj i = . 

e fc , 2 = Re k) , 7fc>2 = Im £*(|i fc) . (28) 



At 0((k F /A) 3 ) we find the equation: 



£fc,3 + «7fc,3 = S 3 (^t, fc) + Ej 



<% 



(29) 



Then the third order contribution to e p and 7p is just given by the sum of real and imaginary 
parts of the diagrams, respectively: 

e M = Rfi^r« 1 (&,*)+IteSS m (^ ) *)+Ite^ 8) (^ > *)+Ite5^ 9) (^ ) A ; ) , (30) 



J 3(6) I 2M ' 'V iLC ^3(7) v 2M ' ^ ^ lbc ^3(8) v 2Af ' 'V ^ iLC ^3(9) v 2M ■ 
"3(6) ( ^) + I m ^3(7) (2M1 



7M = ImE^^ + ImEW^fc) . (31) 



The real and imaginary parts of quasiparticle energy and width in units of the Fermi energy 
Ep = k 2 F /2M up to order {k F /Af are shown in Figs. |3| and [|, respectively. As an example, 
we have chosen the parameters k F a s = k F a p = k F r s = 0.5 and g = 4. The width ^ k shows 
the behavior predicted by the Lehmann representation, that is 7^ > for k < kp and % < 
for k > k F . Our results agree qualitatively with the previous incomplete calculation of Bund 
and Wajntal |1{J where the particle-hole and effective range contributions were omitted. For 
those contributions for which a detailed comparison is possible, quantitative agreement is 
obtained as well. 

Using Eqs. (p8|,30,|31), we can immediately calculate the third order contributions to the 



above mentioned observables. As the effective mass M* is proportional to the inverse of the 
slope of the excitation energy at the Fermi surface, we need the complete result for e& up to 
order (k F /A) 3 : 

W = (l + ( 9 -l)^J 2 (l-7.o g2) + ^ 

+0.11(g-l)(g-3)(k F a s f -0.15(g-l)(k F a s ) 3 ) (32) 
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0, 







0.5 



k/k r 



FIG. 3: The quasiparticle energy e p in units of the Fermi energy Ep up to order (&p/A)° (dotted 
line), (hp/A) 1 (dashed line), (kp/A) 2 (dash-dotted line), and (kp/A) 3 (solid line) for kpa s = 
k F a p = k F r s = 0.5 and g = 4. 



0.15 



0.1 



rjf 0.05 
8^ 



-0.05 




k/k r 



FIG. 4: The quasiparticle width 7 P in units of the Fermi energy Ef up to order (kp/A) 2 (dash- 
dotted line) and up to order (k F /A) 3 (solid line) for k F a s = k F a p = k F r s = 0.5 and 5 = 4. 



In order to check our calculation, we calculate the third order contribution to the chemical 



potential and the energy per particle. According to Eq. (24), the contribution to the chemical 
potential at 0((k F / A) 3 ) is given by: 



M(3) 



kp 
2M 



(g - l)——(k F a s ) 2 k F r s + (g + l)——(k F a 



15tt 



15tt 
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+{g - 1){0.20 + {g- 3)0. 15} (k F a s ) 3 



(33) 



The third order contribution to the energy per particle is then given by: 



-) 



kp r 



2M 



(g - l)—-(k F a s ) 2 k F r s + (g + l)^(k F a p ) 3 



+ (g - 1){0.076 + {g- 3)0.057}(fc F a a 



(34) 



which is in good agreement with previous results |L2|, [L3|, [L4|, [15 . 

Finally, we also calculate the wave function renormalization Z defined in Eq. (|21~D. 3 It 
is related to the discontinuity in the momentum distribution at the Fermi surface f23| . We 
obtain for Z(k F ): 



z(k F )- 1 



7T Z 



log 2 (k F a s 



+ | Q 184 _ _ 3 ) .224}(A; F a s ) 3 



TV 



(35) 



We used the results of Sartor and Mahaux [^(| for the occupation numbers to obtain an 
analytical expression for the second order term. We verified this result numerically. At third 
order, the results were computed numerically. 



V. SUMMARY AND CONCLUSIONS 

In this paper, we have performed the first complete calculation of the the quasiparticle 
excitation energies and widths in the dilute Fermi gas to order (k F /A) 3 . In order to check 
our calculation, we have derived and compared the chemical potential and the energy per 
particle at (k F /A) 3 to previous results. 



The calculation has been carried out using an EFT for dilute Fermi systems [[15|]. Using 
dimensional regularization with minimal subtraction the power counting in this EFT is 
particularly simple. All loop momenta are converted to Fermi momenta which leads to a 



fully perturbative expansion. This is in contrast to previous calculations fL2| , |T3| , |14f . There, 
a more general set of diagrams had to be first summed and then expanded with care to avoid 
double counting. The EFT approach is controlled and allows for systematic improvements 
by simply going to the next order in the power counting. Errors can be estimated reliably 
using dimensional analysis and error plots. 

To order (k F /A) 3 , no three-body input is required. At (9((fcp/A) 4 ), however, a contact 
three-body interaction enters. Without the three-body term, the logarithmic divergence 
appearing at this order can not be renormalized |T3J. A complete calculation of the energy 



density to this order is in progress [25 . 



3 Note that Z depends on the representation of fields and is not an observable in the usual sense. See 
Ref. |24j] for a related discussion on the representation dependence of the momentum distribution. We 
quote Z in the standard representation of fields in which no interaction terms proportional to the classical 
equation of motion are included. This representation is equivalent to the Hamiltonian used in Ref. |20|] . 
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The extension of the present EFT to the nuclear matter problem is not straightforward. 
The breakdown scale is set by the lowest momentum scale in the effective range expansion 
which is of order 10 MeV for nuclear matter because of the large S-wave scattering lengths. 
Similar problems arise for atomic systems near a Feshbach resonance. For scattering in the 
vacuum, the appropriate power counting to deal with the large scattering lengths is known 
ptj| , p7[ . It requires to sum all Cq interactions which for two-particle nonrelativistic scatter- 
ing simply involves summing a geometric series. At finite density, the diagrams involving 
only Cq vertices are not easily summed because both particles and holes are present. To 
make the problem tractable, an additional ordering scheme such as the hole-line expansion 
in the standard nuclear physics approach is needed. Whether the hole-line expansion can be 
justified in EFT using power counting arguments remains to be seen. An interesting possi- 
bility is to exploit the geometry of intersecting Fermi spheres inherent in the finite density 
loop integrals which leads to an expansion in 1/D where D is the space time dimension p8| 



Furthermore, an expansion in the inverse spin-isospin degeneracy g = 4 might prove useful 
||29|| . The nuclear matter problem provides additional challenges. Since hp ~ 270 MeV, the 
pion exchange is long range and cannot be treated within a purely short range EFT. Both 
the large scattering length problem and the inclusion of the long range pion exchange are 
important questions to be addressed by future work. 



Recent progress in the field of cold Fermi gases [16], [L7| also opens the possibility of 
applications of EFT methods in the physics of cold fermionic atoms. In the case of cold 
bosons, EFT methods have already been applied successfully to a number of processes 



31, 32, m 
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APPENDIX A: PROPER SELF-ENERGY INTEGRALS 



We illustrate the calculation of the self-energy diagrams in Section |J| by applying the 
Feynman rules to diagram (7) from Fig. |2|. This diagram has a symmetry factor of 1/4. 
Thus, one obtains the following expression: 

SS ( 7)(P) = H C °) 3 I7^r / d " k J d4 <? / d 4 U 5 GoCk)G (p~q)GoCk + q)G (p-l) . (Al) 
After performing the contour integrations we have: 



-3(7) 



(P) 



cf 9 ~ l 



d*k / d 3 g / dH 



(2rr) 9 

20(fc - fc F )fl(|k + 1| - fc F )0(|p - 1| - £; F )0(fc F - |k + q|)0(fc F - |p - q|) 

(p + uj k - ^k+ q - ^ P -q - ze)(u; k+q - u; k+ i - u; p _i + cj p _ q + ie) 
6(k - k F )6{k F - |k + l|)0(ifc F - |p - l|)0(Jfe F - |k+ q|)0(fc F - |p - q|) 



p-q 



ie 



(A2) 



20(£; F - fc)0(|k + 1| - k F )6(\p - 1| - A; F )0(£; F - |k + q\)6(k F - |p - q|) 



(p + - ^ k +i - Wp_i + ie) (cu k+ i - ^ k+q + ^p_i - uj- 



p q 



it 



6{k F - fc)fl(|k + 1| - k F )9(\p - 1| - fc F )g(|k + q| - k F )6(\p - q| - fc F ) 
(p + w fc - w k+ i - u; p _i + ie) (p + u k - w k+q - Wp_q + ie) 

By setting this expression on-shell (that is po = p 2 /(2M)) and an appropriate substitution 
to dimensionless variables using 



k = k F u, q = k F -(y — u — 2t), 1 = fc F -(v — u — 2r), p = A; F v 



(A3) 



Eq. 



can be simplified significantly: 



J 3(7) 



(P) 



c 3 M 2 fc F (p - 1; 

(2tt) 9 
~20(u- l)6»([s - 



d 3 M / dH / d 3 r 



l)0(|s 



1)0(1 



t|)0(l 



t|) 



(2q 2 - 2t 2 - ie)(2t 2 - 2r 2 - ie) 
6(u-l)6{l - |s-r|)fl(l - |s-r|)fl(l - |s — t|)6>(l - |s + t|) 

{2q 2 - 2r 2 - ie){2q 2 - 2t 2 - ie) 
20(1 - u)6{\s - r| - l)6>(|s - r| - 1)0(1 - |s - t\)6{l - |s + t|) 
{2q 2 - 2r 2 - ie){2r 2 - 2t 2 - ie) 
1 - u)9(\s - r| - l)6{\s - r| - l)0(|s - t| - l)0(|s + t| - 1) 



(A4) 



+ {2q 2 -2r 2 + ie){2q 2 -2t 2 

We consider only the first term in the square brackets: 



ie) 



d 6 u / d 6 t / d 6 r 



3 29(u- l)0(|s-r| - l)0(|s-r| - 1)0(1 - |s-t|)0(l - |s + t|) 



(A5) 



(2q 2 - 2t 2 - ie)(2t 2 - 2r 2 - ie) 

where the constant overall prefactor in Eq. (|A4|) has been dropped. By analyzing the step 
function, we can rewrite the integral as 



It = 8tT 



dx 



du I dz 
o 



z- (s,z) 



dt / dy 



dr 



u 2 t 2 r 2 e(l - s) 



2 +( s >: 



(q 2 - t 2 - ie)(t 2 - r 2 + ie) 



(A6) 
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with 

z±(s, z) = ±SZ + yjl- s 2 (l- z 2 ) . (A7) 
The integral over r contains a linearly divergent term which can be separated by writing 

poo ^2 rz+(s,y) ^ 

dr- = / dr- + / ' dr- . (A8) 

z+(s, y ) t 2 — r 2 + ie Jo t 2 — r 2 + ie Jo t 2 — r 2 + ze 



The ultraviolet divergence is contained in the first term on the right-hand side of Eq. 
In dimensional regularization with minimal subtraction the linear divergence is subtracted 
automatically and the finite part of this term is purely imaginary. Thus, we can write X\ as: 

if 1 , f°° , f 1 , r- (s,z) , u 2 t 2 e(i - s) ( mt\ , K , 

-s^ f a. r d« r a, / M *" ) <u /' d V r +< " y) dr- " 2f2r2s(i - s) 



i 



(g 2 -t 2 -ie)(t 2 -r 2 + ze) ' 



Separating the real and imaginary part of this integral and some further straightforward 
manipulations lead to expressions that can be integrated numerically The remaining terms 
in the square brackets in Eq. (|A4| ) are treated analogously. 
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